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I. INTRODUCTION 

Jordan-Brans-Dicke (JBD) theory 0,0] is the simplest 
extended theory of gravity, depending on one additional 
parameter, the Brans-Dicke coupling u), as compared to 
General Relativity. As Einstein's theory is recovered in 
the limit uj —> 00, there will always be viable JBD theories 
as long as General Relativity remains so too. As such, it 
acts as a laboratory for quantifying how accurately the 
predictions of General Relativity stand up against ob- 
servational tests. The most stringent limits are derived 
from radar timing experiments within our Solar System, 
with measurements using the Cassini probe [j| now giv- 
ing a two-sigma lower limit uj > 40, 000 (improving pre- 
existing limits by an order of magnitude). 

With precision cosmological data now available, par- 
ticularly on cosmic microwave background (CMB) 
anisotropics from the Wilkinson Microwave Anisotropy 
Probe (WMAP) @, it has become feasible to obtain com- 
plementary constraints from the effect of modified grav- 
ity on the structure formation process, as suggested in 
Ref. 0. That paper focussed on the way that u alters 
the Hubble scale at matter-radiation equality, which is 
a scale imprinted on the matter power spectrum, in an 
attempt to identify how large an effect can be expected. 
Subsequently, the expected total intensity and polariza- 
tion microwave anisotropy spectra in the JBD theory 
were computed, and a forecast of the sensitivity to uj of 
data from the WMAP and Planck satellites carried out 
exploiting a Fisher matrix approach |7|. 

In this paper we make a comprehensive comparison 
of predictions of the JBD theory to current observa- 
tional data, using WMAP and other CMB data plus the 
galaxy power spectrum as measured by the two-degree 
field (2dF) galaxy redshift survey. We define JBD models 
in terms of eight parameters, which are allowed to vary 
simultaneously. Our paper is closest in spirit to work by 
Nagata et al. [8j, who considered a more general model, 
the harmonic attractor model, which includes JBD as a 
special case. However their dataset compilation was re- 



stricted to the WMAP temperature power spectrum. 

The constraint we will obtain is not competitive 
with the very stringent solar system bound given above 
(though the analysis of Ref. [Tj indicates that a limit as 
high as 3000 might eventually be reached by the measure- 
ments of the Planck satellite), but it is complementary in 
that it applies on a completely different length and time 
scale. Such constraints can therefore still be of interest in 
general scalar-tensor theories where uj is allowed to vary; 
for instance Nagata et al. |8( find that in some parameter 
regimes of the harmonic attractor model the cosmologi- 
cal constraint is stronger than the Solar System one. In 
that regard, our result is most comparable to cosmologi- 
cal constraints imposed on uj from nucleosynthesis, which 
give only a weak lower limit of uj > 32 [|| . 

II. FORMALISM 
A. Background cosmology 

The Lagrangian for the JBD theory is 

2 

where the Brans-Dicke coupling uj is a constant, and Q(t) 
is the Brans-Dicke (BD) field whose present value must 
give the observed gravitational coupling. We have in- 
cluded factors of mp\ to define $ as dimensionless. 

The e quat ions for a spatially-flat Friedmann universe 
are dHia 

f a\ 2 a $ uj f$\ 8tt 

UJ = 6 {*) + 3mf 1 ¥ P; (2) 

* + 3 ~ * = to TT\ 2 (P " 3P) > ( 3 ) 
a (2uj + 3)mpj 

where a(t) is the cosmological scale factor, and p and 
p are the energy density and pressure summed over all 
types of material in the Universe. 
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FIG. 1: Evolution of the BD field from early in radiation 
domination to the present. It is just possible to see the evolu- 
tion of $ increase as A domination sets in. The cosmological 
parameters are ui = 200, Ho = 72, and p mi o = 0.3 in units of 
the standard cosmology critical density. 



The Universe is assumed to contain the same ingredi- 
ents as the WMAP concordance model [Bj], namely dark 
energy, dark matter, baryons, photons and neutrinos. We 
make the simplifying assumptions of spatial flatness, dark 
energy in the form of a pure cosmological constant, and 
effectively massless neutrinos whose density is related to 
that of photons by the usual thermal argument. The 
present value of $ must correctly reproduce the strength 
of gravity seen in Cavendish-like experiments, which re- 
quires 
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where here and throughout a subscript '0' indicates 
present value. We will assume that the value of $o m 
our Solar System is representative of the Universe as a 
whole, though this may not be absolutely accurate [Tl| . 
We also assume that the initial perturbations are given 
by a power-law adiabatic perturbation spectrum. 

When the Universe is dominated by a single fluid there 
are a variety of analytic solutions known , where $ is 
typically constant during a radiation era, slowly increas- 
ing during a matter era, and then more swiftly evolving 
as dark energy domination sets in. However we need 
solutions spanning all three eras and so will solve the 
equations numerically, for which we use the integration 
variable N = lnet/ao. An example of the evolution is 
shown in Fig. ^ 

The basic parameter set we use to build our cosmolog- 
ical models contains the following parameters 



u> Brans-Dicke coupling 

Hq present Hubble parameter [kms _1 Mpc 

Pb baryon density 

pc cold dark matter density 



curvature perturbation amplitude 
ng perturbation spectral index 
t reionization optical depth 
b galaxy bias parameter, P g g/P mm 

where b 2 = Pg g /P m m is the ratio of the (observed) 
galaxy power spectrum to the (calculated) matter power 
spectrum. Other parameters are fixed by the assump- 
tions above, and the radiation energy density is taken as 
fixed by the direct observation of the CMB temperature 
T = 2.725K 0. 

An important subtlety that must be taken into account 
is that the extra terms in Eq. (2), plus the Cavendish- like 
correction to the present value of <&, means that the usual 
relation between the Hubble parameter and density, used 
to define the critical density and hence density parame- 
ters, no longer applies. Generically, the extra terms re- 
quire an increase in the present value of p to give the 
same expansion rate, the correction being of order 1/lo. 
Because of this subtlety, we define the density parameters 
^b,c by dividing by the critical density for the standard 
cosmology, meaning that the density parameters don't 
quite sum to one for a spatially-flat model. 

Operationally, we proceed as follows. We seek a back- 
ground evolution corresponding to a particular value of 
h = i/o/100 and of the present physical matter density. 
We can assume the initial velocity of the BD field <£> is 
zero deep in the radiation era, which leaves us two pa- 
rameters, the early time value of <3? and the value of the 
cosmological constant, to adjust in order to achieve the 
required values. This is a uniquely-defined problem, with 
the necessary values readily found via an iterative shoot- 
ing method. 



B. Perturbation evolution 

We carry out the evolution of density perturbations us- 
ing a modified version of the code DEFAST, based on CMB- 
fast |3] and originally written to study quintessence 
scenarios where the dark energy scalar field is minimally 
[TH or non-minimally 01 coupled to the Ricci scalar. 
The architecture of DEFAST is based on the version 4.0 
of CMBFAST, although there has been a progressive code 
fork in the subsequent versions. DEFAST takes as input 
the parameter set described in the previous subsection, 
and returns the microwave anisotropy spectra (for tem- 
perature and polarization) and the matter power spec- 
trum. A dynamical and fluctuating scalar field, playing 
the role of the dark energy and/or the BD field, is in- 
cluded into the analysis together with the other cosmo- 
logical components, following the existing general scheme 



In order to bring the model description into the for- 
malism used by DEFAST, we redefine the BD field and 
coupling according to 
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which brings the Lagrangian into the form 

C = - ~9„<OT + Anatter , (6) 

where <j) is now a canonical scalar field non-minimally 
coupled to gravity. We implement the cosmological con- 
stant in the code by giving a constant potential energy. 

Our calculations include the effect of perturbations, 
with the initial perturbations in </> fixed by the require- 
ment of adiabaticity. The correction to the background 
expansion rate from the dynamics of </> is the most rel- 
evant effect on the CMB power spectrum, appearing as 
a projection plus a correction to the Integrated Sachs- 
Wolfe (ISW) effect, as discussed in detail in Ref. [l6|. 



C. Data analysis 

The data we use are taken from WMAP [2l| and the 
2dF galaxy redshift survey expressed as 32 bandpowers 
in the range 0.02 < k < 0.15/i _1 Mpc H^. In order to 
incorporate the 2dF data, the galaxy bias parameter b 
is taken to be a free parameter for which the analytic 
marginalization scheme of Ref. [23 can be applied. We 
also consider the effect of including the high-^ CMB data 
from VSA [3, CBI [3, ACBAR UJ. 

Our present analysis does not include supernovae data. 
Inclusion of the modification to the luminosity distance 
from uj would be straightforward. However the variation 
of the gravitational coupling G means that supernovae 
can no longer be assumed to be standard candles, and 
Ref. [24| suggests that the effect from varying G domi- 
nates. Further, inclusion of supernovae data may be par- 
ticularly susceptible to the possibility that the local value 
of <E> in the vicinity of the supernova may not match the 
global cosmological value [Tl| ■ Nevertheless, it would be 
interesting to investigate robust methods for including 
such data, also in connection with alternative observa- 
tional strategies p5| . 

We carry out the data analysis using the now-standard 
Markov Chain Monte Carlo posterior sampling tech- 
nique, by modifying the June 2004 version of the Cos- 
moMC program |2f| to call DEFAST to obtain the spec- 
tra. CosmoMC computes the likelihood of the returned 
model and assembles a set of samples from the poste- 
rior distribution. We take full advantage of CosmoMC's 
MPI capabilities by running the code across 19 Sun V60x 
Xeon 2.8GHz processors. The Metropolis-Hastings algo- 
rithm is run at a temperature of 1.3 in order to better 
sample the non-Gaussian direction of our posterior dis- 
tribution which results from the degeneracy between Hq 
and ln£, both of which have a strong effect on the angu- 
lar diameter distance. The final chains are then cooled 
and importance sampled |2(j. It can be noted that for 
the purposes of posterior sampling, we have parametrized 
the JBD cosmology using In £ = — In Auj simply because 
it is more straightforward to obtain the samples we need, 
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FIG. 2: Marginalized ID posterior distributions (solid lines) 
on the base parameters as listed in Section II. Also displayed 
is the mean likelihood of the binned posterior samples (dotted 
lines). 



while simultaneously suppressing the possibility of jump- 
ing to regions with uj < 1. Specifically, we use a flat 
prior on ln£ 6 [—9,-3] where the lower cutoff has been 
adjusted to the point where the likelihood function is no 
longer sensitive to the effect of varying the Brans-Dicke 
parameter and the ACDM model is thereby recovered. 
As usual, this Jeffreys prior, which is defined here as a 
flat prior on the logarithm of a parameter of unknown 
scale, has the interestin g p roperty of invariance under 
scale reparametrizations 27] . For this reason it serves as 
a reasonable substitute for working with a more desir- 
able physical parameter which could be identified to iso- 
late and give a linear response in the ISW effect, mainly 
responsible for the upper bound on ln£. 

The optical depth r is parametrized using Z = 
exp [— 2t], where Z 1 ' 2 is the fraction of photons that re- 
main unscattered through reionization, since the combi- 
nation Ag Z is well constrained by the CMB. 

The results that we present are based on around 
100,000 raw posterior samples, and while the basic con- 
straints can be derived with significantly fewer samples, 
this large number assures more robust constraints on the 
derived parameter to when we use importance sampling 
in order to adjust for the change in prior density |2q. 



III. OBSERVATIONAL CONSTRAINTS 

Turning first to the constraints on the basic param- 
eter set, from Figure |21 we note the overall consistency 
of our results with the current observational picture (see 
for example Ref. [5( and a work by two of the current au- 
thors Ref. 28]), finding the 99% marginalized probability 
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FIG. 3: Marginalized ID posterior distributions (solid lines) 
on the BD parameter ln£ (upper panel). Also displayed are 
the derived importance sampled constraints (correcting for 
the change in prior density) on the more familiar lu (lower 
panel, no smoothing). We obtain a 95% marginalized proba- 
bility bound of ln£ > —6.2, corresponding to a bound on the 
BD parameter lu > 120. 

regions to be 

0.021 < n B h 2 < 0.027, 0.10 < fl c h 2 < 0.15 , 

61 < H < 80, 0.57 < Z < 0.97 , (7) 
0.92 < n s < 1.07, 19 < A s < 33. 

Note that part of our allowed region lies outside the priors 
assumed by Nagata et al. 8]. As usual for joint analyses 
of CMB and galaxy power spectrum data, it is unneces- 
sary to impose a further constraint on Hq. 

The primary focus of our study has been to derive con- 
straints on the BD parameter for which, from the outset, 
we have expected only to find a one-sided bound; the sit- 
uation can only become more interesting when both the 
angular diameter distance and the recombination history 
become much better probed by the CMB. This expecta- 
tion is indeed confirmed by the data, as shown in Figure|3| 
in which we display the region of highest posterior den- 
sity. The lower panel detailing the posterior constraint on 
to has been obtained by importance sampling to correct 
for the change in prior density when changing param- 




FIG. 4: Marginalized 2D posterior distribution in the ln£-i/o 
plane. The solid lines enclose 95% and 99% of the probabil- 
ity. Under this parametrization there is clearly a geometrical 
degeneracy. 

eters from ln£ to lu (we note that the mean likelihood 
of the binned posterior obtained from sampling ln£ per- 
forms well for putting a bound on u, demonstrating less 
sensitivity to the details of the prior density). 

We obtain calculate marginalized probability upper 
bound and the main result of this paper to be 

ln£ < -6.2, 95%, 

ln£ < -5.7, 99%. (8) 

The corresponding marginalized probability lower 
bounds on the BD parameter are found to be 

lu > 120, 95%, 

lu > 80, 99%. (9) 

This bound is nicely consistent with the expectation for 
WM AP given by the Fisher matrix analysis of Ref . • 

We present in Figure 0] the 2D posterior constraints 
in the ln^-_ffo plane, in order to demonstrate the 
degeneracy and covariance between these two param- 
eters. In a more refined analysis, one could re- 
place Hq with the dimensionless parameter r s /DA 
more appropriate to the study of the CMB, where 
r s is the sound horizon at recombination and Da is 
the angular diameter distance to the last-scattering 
surface |2^. Finally, in Figure we display two 
models, our best-fit ACDM model with parameters 
0={n B h 2 , n c h 2 ,H , Z, ns, 10 1Q A S , w}= {0.023, 0.12, 66, 
0.79,0.96,23.2,oo}, and a best-fit JBD model with pa- 
rameters 6= {0.024,0.13,79,0.80,1.03,24,70}, in order 
to illustrate how the observables change at finite lu. Here 
the JBD model lies in the vicinity the contour enclos- 
ing 99% of the posterior probability distribution and was 
selected by running a short Monte Carlo exploration at 
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FIG. 5: A comparison between a ACDM model (solid line) 
and a JBD ACDM model with uj = 70 (dashed line). The 
data are the 2dF galaxy power spectrum and the models the 
matter power spectrum convolved with the 2dF window func- 
tions, and whose overall amplitude is left as a free parameter. 
Detailed parameters are given in Section ITTT1 

fixed uj = 70. Note that although in principle the param- 
eter ln£ could be extended to — oo, whereby the bulk of 
the parameter space would be composed of the ACDM 
model, in practice it is reasonable to adjust the lower cut- 
off to the point where the likelihood function loses sen- 
sitivity to the variation of ln£ so that the Brans-Dicke 
model alone is explored by the MCMC. Consequently, 
the probability contours can reasonably be interpreted 
to describe the most credible region of the Brans-Dicke 
model parameter space. 

Our current analysis leaves the bias parameter free, 
and so constrains only the shape of the matter power 
spectrum. We note however that the JBD model has a 
significantly higher amplitude, indeed requiring a mod- 
est antibias b ~ 0.98, which at least in part is due to the 
more rapid perturbation growth (5 oc a 1+1 ^ during mat- 
ter domination 0) in the JBD theory. For comparison 
the ACDM model has a best-fit bias b = 1.2. This sug- 
gests that precision measures of the present-day matter 
spectrum amplitude, as for instance may become avail- 
able via gravitational lensing, could significantly tighten 



constraints. We also note that there is a shift in the 
location of the baryon oscillations in the matter power 
spectrum as compared to the ACDM model; these are 
mostly erased by the 2dF window function, 1 but future 
high-precision measurements of those may also assist in 
constraining uj. 

We have carried out the same analysis including also 
the data from VSA, CBI and ACBAR in the multipole 
range 600 < I < 2000. This high-£ data leads to a slightly 
tighter bound on the Brans-Dicke parameter, ln£ < —6.4 
corresponding tow > 177 at 95% marginalized probabil- 
ity. However, at the same time inclusion of this new 
data leads to an unexpectedly large shift in the spec- 
tral index, to 0.90 < ng < 1.00 at 95% marginalized 
probability, so that the Harrison-Zel'dovich spectrum is 
only just included (this statement remains true in the 
general relativity limit). Whether this points to some 
emerging tension in the combined dataset, a harmless 
statistical fluctuation, or a hint of the breaking of scale- 
invariance, can be addressed only in the light of the next 
round of CMB observations. While our constraint on ln£ 
marginalizes over ns , in the interests of quoting a robust 
bound we have given as our main result the weaker limit 
obtained without including the high-^ data. 

Our ultimate constraint uj > 120 can be compared with 
that of Nagata et al. who quote results correspond- 
ing to uj > 1000 at two-sigma and uj > 50 at four-sigma. 
The former constraint is much stronger than projected in 
Ref . , and stronger than one would expect from a naive 
assessment that the corrections to observables should be 
of order 1/uj, If we plotted a model with uj — 1000 in our 
Figure it would lie practically on top of the ACDM 
model. However their latter constraint is in reasonable 
agreement with ours, and they do highlight that it is this 
constraint which corresponds to a sharp ridge of deterio- 
rating chi-squared in their analysis, indicating that their 
constraint should conservatively be taken as uj > 50. 



IV. CONCLUSIONS 

We have derived a constraint on Jordan-Brans-Dicke 
gravity from current cosmological observations, including 
cosmic microwave background (CMB) anisotropy data 
and the galaxy power spectrum data. Our main result 
is to obtain a 95% marginalized probability lower bound 
on the Brans-Dicke parameter uj > 120. This result is 
complementary to the very strong Solar System limit 
provided by Cassini, uj > 40000, as it probes entirely 
different length and timescales. Our analysis is based on 



Our analysis used 2dF data from Percival et al. |22| . preced- 
ing the more recent 2dF data analysis which shows evidence of 
baryon oscillations We would not expect inclusion of this 

new data (not yet publicly available) to significantly change our 
results. 
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a Markov Chain Monte Carlo technique varying the basic 
cosmological parameters and 10. 

At the present precision level, the greatest part of the 
constraining power comes from the shape of the CMB 
acoustic peaks, in particular from the first-year obser- 
vations of WMAP. Therefore, assuming an extension to 
four years of the WMAP observations, we expect some 
further improvement on the limit on u> from cosmology. 
Further help is also expected from other structure for- 
mation data, as they improve quality and precision in 
coming years. In particular we have highlighted that an 
accurate measure of the present-day matter power spec- 
trum amplitude, for instance from gravitational lensing, 
may be powerfully constraining when compared to the 
primordial amplitude from the CMB. 

A leap forward in this and other contexts is expected 
from the observations of the Planck Surveyor probe, to 
be launched in 2007. Those observations are expected 
to be cosmic variance limited for the whole spectrum of 
CMB temperature anisotropy down to the damping tail, 
and to provide an accurate measurement of the gradi- 
ent mode of the CMB polarization and its correlation 
with total intensity up to the sixth acoustic peak |3l| . 
According to the forecasts of Chen and Kamionkowski 



, the limit on u> from Planck should be around an or- 
der of magnitude stronger than that from WMAP, and 
hence vastly stronger than the nucleosynthesis constraint. 
Whether that improvement can be realised from actual 
Planck data, of course, remains to be seen. 
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